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Annotation: Ensuring the strength, reliability and efficiency of hydraulic structures requires 
comprehensive consideration of the static and dynamic loads acting on the structure in a single 
statement. When calculating the seismic effects of structures such as ground dams that perceive 
water pressure, in addition to seismic inertial forces, it is necessary to take into account the 
hydrostatic and hydrodynamic pressure of water on the pressure face, as well as the pore pressure 
in the body of the dam as a multiphase medium. Currently, consideration of seismic water pressure 
is mandatory when designing hydraulic structures, because our republic is located in a seismically 
active area. 


The purpose of this article is to describe a methodology for taking into account the hydrodynamic 
pressure of a liquid during seismic impacts on a dam without fastening plates, taking into account 
the interaction of the liquid and the skeleton of the soil. 


Keywords: earth dam, hydrodynamic pressure, dynamic impact, stability, oscillation period, 
seismic active area, water-supporting hydraulic structure. 


Introduction. The existing system of designing and monitoring the state of hydraulic structures has 
significant shortcomings, consisting in the absence of reliable and objective criteria that 
characterize the safety of structures. Therefore, when designing, it is necessary to choose the type 
of structure, the predominant oscillation period of which differs most from the oscillation period of 
the base [1,2]. rrently, consideration of seismic water pressure is mandatory when designing 
hydraulic structures, because our republic is located in a seismically active area. The magnitude and 
distribution of the seismic water pressure acting on the pressure face of the structure, as well as the 
pore pressure inside the body of the dam depends on many factors: the direction, intensity and 
regularity of the seismic impact, the geometric configuration of the structure and the reservoir, the 
elastic compliance of the structure and the base, the ratio of the prevailing frequencies of seismic 
impacts and the natural frequencies of the structure and the liquid (as acoustic environment), etc. 
[3,4]. The degree of accuracy with which these factors should be taken into account depends on 
what proportion of the total complex of loads acting on the structure is seismic water pressure. 
When determining the hydrodynamic pressure of a liquid, it is necessary to take into account the 
type of construction of a soil dam. This formulation of the problem will make it possible to more 
realistically describe the stress-strain state of ground dams [9, 10]. 


The state of the problem. In accordance with the requirements of KMK 2.06.05-98 “Dams from 
soil materials”, when designing earth dams of classes 1 and 2, justifying calculations of the 
structure structure for basic and special combinations of loads during construction and operation 
should be performed, including: 
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> stress-strain state of the dam; 
> stability of the dam. 


Calculations of all hydraulic structures, foundations and coastal slopes both in the alignment of the 
structure and in the reservoir area should be carried out on static loads determined by linear spectral 
theory (LST) under the influences set by the conditional accelerations of the base A (expressed in 
fractions of g — acceleration of free fall). 


Calculations must be made in accordance with the requirements of the heads of the KMK for the 
design of hydraulic structures of certain types. The calculations should take into account seismic 
loads from the mass of the structure, the attached mass of water (or hydrodynamic pressure), from 
waves in the reservoir caused by an earthquake, and from the dynamic pressure of the soil [5]. 


For water-retaining hydraulic structures of class I, when they are located in areas of seismicity J>7 
points, it is also necessary to make calculations for seismic effects specified by accelerograms 
(velocigrams, seismograms) of earthquakes. 


In this case, the designing organization develops special technical conditions in which: 


> the scheme of setting the seismic impact, the relationship between accelerations in different 
directions, changes in amplitudes along the perimeter of the canyon, the time shift of the impact 
at various points of the base of the structure are determined; 


> the main provisions of computational or experimental studies are established; 


> the composition and methodology of experimental field and laboratory studies of the physical 
and mechanical properties of the soils of foundations and building materials of structures are 
determined within the framework of the accepted models; 


> criteria are formulated to determine the maximum bearing capacity of the structure and to 
assess its seismic resistance [6]. 


The solution to the problem. A typical reservoir area is shown in Fig.l, and its contour is 
determined by the water level in the reservoir, the outline of the upper face of the dam and the 
accepted length of the calculated area. With a sufficiently elongated area, the boundary conditions 
along the AB line do not significantly affect the results of the joint calculation of the reservoir- 
structure system. 


2 


Fig.1. Finite difference grid for calculating the wave equation in the reservoir: 


points for determining velocities in the dam; o - points for determining stresses in the dam; 
X - points for calculating the potential @ in the water storage. 
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Let the fluid motion equation have the form 
0°D A CD _ p a@ 
Gx? ay? M, Gt? 


(1) 


where ( - is the density of water; Mp, - modulus of elasticity of water; ® - potential speeds: 
“= d@/ dx;v =d®/dy; u,v - vertical and horizontal fluid velocities. 


The hydrodynamic pressure P is related to the velocity potential by the ratio 


O@ 
P=p— (2) 

Ot 
To solve equation (1) along the contour of the computational domain, either the values of the 
potential or the velocity along the normal to the boundary must be specified. The initial conditions 


will be considered homogeneous, that is, for f =0, P(x, y,0)=0 u d®/ dt(x,y,0)=0. 


For simplicity, let's consider the case when O@/ Oy =O (symmetry conditions) are set on AB. 


Wave formation on the free surface will be neglected. The error of the assumption decreases with 
the increase in the height of the dams and is negligible for a dam with a height of /00-200 m. It 
should be noted that at water flow rates characteristic of earth dams, the additional hydrodynamic 
pressure is very small and therefore it does not make sense to determine it with high accuracy. We 
will consider the section of the BC contour to be fixed and from the condition that there is an 
P=O on it, we will set the boundary condition of the @cg =O. Along the AFEDC line, the 


speeds of water movement along the normal to the contour are given. According to the base AF, 
these speeds are determined by the given law of displacement of the base, and along the upper face 
of the dam FEDC, the speeds of soil movement, determined by the joint parallel calculation of the 
"dam - reservoir" system. 


Let us take the structure of the difference grid shown in Fig.l. The physical meaning of equation 
(1) is that the amount of water flowing out of a unit of volume per unit of time (the left side of the 
equation) is proportional to the decrease in pressure due to elastic deformations of water. Consider 
the volume element shaded in Fig. 2, a, with the area: 


O, = (x, +X; KY: a a 14. 


The velocities at the points lying in the middle of the sides will be determined by the formulas: 


. Pi 445 “O74 . 

Ue >, a= — a — 4 
Had Xx; 
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Then the amount of liquid flowing from the volume element in question in | s is equal to: 


2 
ses? ee ; i p O®, 


O..=| u —u ——__—— +] -¥~ = —— (4) 
U fe es _i, res a 2 

i) pd 2 bits BI 2 M, Ot 
Equation (4), when (3) is substituted into it, is a difference analog of equation (1) and can be 
obtained from it by a formal difference approximation of the second derivatives of ® with respect 
to spatial coordinates. However, when considering points on the contour of a region, it is more 
convenient to keep in mind this physical meaning of difference equations. 


Consider the volume element adjacent to the slope of the dam (Fig. 2, b), with the area: 


Oo, = ee PO dl 


¢) 


Fig.2. Difference analogues of wave equations: a - internal point of the reservoir area; b - point on 
the sloping boundary; c - point on the conjugation of the horizontal and inclined sections of the 
contour. 


Let the speed u* be given along the slope along the normal to the boundary. The volume of water 
flowing out of the element, 


aoe a XxX, +X 


where 
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re ae 
. = ee 
j 2sina X. +x. 


YAN 


Let the speed u' given on the contour be represented by its projections “1” and vy": 


«k os ae “0 
uo =—-u SINaA+V COSA 


Then the amount of water considered by the point we get: 


VUE e's X +X. 

ues =? J J 4+) 1 = 

J 2 Z 
Finally , for the point under consideration, we get: 
2 
gle Ge x +x. : OP. 
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2 


As an example illustrating the obtaining of difference analogs of equation (1) at contour points, 
consider the case of conjugation of a horizontal and an inclined section of the contour (points of 


J 
type F and D in Fig.1). The considered element of the area with area @ = 16 “{ Vory. | 
ij i i J7- 


is shown in Fig.2.2, b. 
For this case 


0°@D 


ee ae x p : 
Q =|i ui" Et fay \p-0 —_—__! (6) 
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Similarly, difference equations are compiled for horizontal sections of the contour (the base, the 
berm of the dam) and for points of type E. Thus, the general record of difference equations has the 
following form: 


7D 
=o 4 (7) 
ij i M , ar? 
where 0. - is the flow rate of the fluid flowing from the element in question; O,; - is the area of 
this element. 


Writing in finite differences 


Published under an exclusive license by open access journals under Volume: 3 Issue: 2 in Feb-2023 
Copyright (c) 2023 Author (s). This is an open-access article distributed under the terms of Creative Commons 
Attribution License (CC BY).To view a copy of this license, visit https: //creativecommons.org/licenses/by/4.0/ 


106 


EEE SSS ee eee ee ee 
Journal of Ethics and Diversity in 
JEDIC 


International Communication 


| e-ISSN: 2792-4017 | www.openaccessjournals.eu | Volume: 3 Issue: 2 


O-@ p't4t_os@! 4 @'~4t 
ij ij ij ij 


and expressing explicitly the value of oe , we get the recurrent formula: 


M 
p+ 4-20! —@'- 44 At? ame 
q] Y Y] p (8) 


allowing step by step to calculate the values of @; - (and hence P) for any t>0. The joint 
calculation of the dam-reservoir system is carried out according to the following scheme: 


1) let at some time step ¢ be known: 


a) the value of @’ and @' in the area occupied by the reservoir; 
b) the values of the speeds of soil movement in the dam; 


2) according to the set speeds on the contour and (8) in the entire region, the value @ is 


determined at the next time Dea : 


3) along the contact line of the dam with the reservoir, the pressure is determined by the formula 
op 4t_@! 


Pag 9) 
Z=P Ai ( 


which is the difference analogue of equation (2); 


4) pressure P; is used in the formulation of equations of dynamic equilibrium in the contour 


elements of the dam, which determine the speed of the soil at time f+ At. The process is then 
repeated for the next time step. 


The speed of water movement cku and ckv at the contact of the dam with the reservoir is 
determined by the formulas 


a? =u, + Wm; - 
Ded che (10) 
vo =v, + Wym; 


where ™ - is porosity; 4“ ,v - components of the skeletal material velocity vector; W ,w_ - 
c oc xy 


components of fluid velocities relative to skeletal material. 


Conclusion. The developed calculation method assumes the determination of the stress-strain state 
of the ground structure at any stage of deformation and therefore primarily includes an assessment 
of the work on the permissible deformations of the elements of the structure. Local destruction at 
some point of the ground does not mean that the structure is in a dangerous position. Only in the 
case when the whole area, which has an outlet to a free surface, goes into a limiting state, plastic 
deformations in this area will accumulate without a boundary. The progressive nature of the 
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accumulation of plastic deformations in any element of the structure indicates that the structure is 
in a state of extreme strength. 


Thus, the criterion for the overall stability of a soil structure is the absence in this structure of areas 
with a progressive accumulation of plastic deformations. In the numerical analysis of the stress- 
strain state of the designed dams, we used those introduced in the definitions of the limit state of 
the soil structure. 
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